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Abstract. We use the results of our recent NLTE model 
atmosphere analysis of central stars of old planetary neb- 
ulae (PN) to calculate distances. We perform a compar- 
ison with three other methods (trigonometric parallaxes, 
interstellar NaD lines, and Shklovsky distances) and dis- 
cuss the problem of the PNe distance scale. The result 
of the comparison of our spectroscopic distances with the 
trigonometric distances is that the spectroscopic distances 
are 55% larger. Since using trigonometric parallaxes with 
large relative measurement errors can introduce system- 
atic errors, we carried out a Monte Carlo simulation of the 
biases introduced by selection effects and measurement er- 
rors. It turns out that a difference between both distance 
scales of the observed size is expected for the present day 
data if the underlying distance scales are identical. Thus 
our finding is essentially a confirmation of the spectro- 
scopic distance scale! Good agreement is found between 
the spectroscopic distances and distances derived from the 
interstellar NaD lines. 

All three independent methods of distance measure- 
ment indicate that the widely used "statistical" distance 
scales of the Shklovsky type are too short for old PNe. 
A correlation with nebular radii exists. The most likely 
explanation is an underestimate of the nebula masses for 
large PN. Implications for the nebula masses are discussed. 
Estimates of the PNe space density and birthrate, which 
are based on Shklovsky type distances, therefore give too 
large values. 

Key words: stars: distances — planetary nebulae: gen- 
eral — white dwarfs — stars: fundamental parameters 



1. Introduction 

Distances of PNe are usually difficult to determine and 
are a long standing problem of PNe investigations. The 
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vast majority o f publ ished distances of PNe are based on 
the Shklovsky (|l956|) method or derivatives of it (often 
called statistical distances). This method allows the cal- 
culation of PN distances from the measurement of the re- 
combination line H/3 and the angular diameter. However, 
Shklovsky distances are notoriously smaller than distances 
derived from model atmosphere analysis of the central 
stars (cf. Mendez |l988l ). 

The question of PNe space densities and birth rates 
is closely rel ated t o the distance scale problem. Ishida & 
Weinberger ( 1987 ) compiled a list of nearby PNe, which 
contains mostly old, evolved nebulae (actually many of 
our PNe were selected from this list). Ishida & Weinberger 
collected distance determinations from literature and com- 
puted the space density and birth rates of this local sample 
of PNe. The derived birth rate of 8 ■ 10~ 12 yr _1 pc~ 3 is too 
high to be in accordance with estimates of white dwarf 
birth rates (2.3 • 10~ 12 yr _1 pc~ 3 ; Weidemann 1991). Since 
every central star of a PN (CSPN) should become a white 
dwarf this yielded a real dilemma. Taken at face value this 
would indicate that current white dwarf samples are very 
incomplete and the white dwarf birth rates are seriously 
underestimated. A certain fraction of white dwarfs may be 
hidden in binaries, indeed. Weidemann ( |1991| ) used a very 
local sample of white dwarfs (d < 10 pc) for his estimate 
and applied corrections for incompleteness and binarity. 



Pottasch (1996) reevaluated the PN space density and de- 
rived a value of 3 • 10~ 12 yr^pc -3 , lower than Ishida & 
Weinberger's result, but still higher than the estimated 
white dwarf birthrate. The difference in the PNe birthrate 
can be traced back to the use of different distance scales. 
While Ishida & Weinberger's collection is mainly based on 
statistical distance estimates (mainly Shklovsky distances 
and derivates of this method), Pottasch excluded statisti- 
cal distance determinations. 



In Paper IV of this series (Napiwotzki 1999) we pre- 



sented the results of an NLTE model atmosphere analysis 
of 27 central stars of old PNe. This analysis of a reason- 
ably sized sample of central stars of old PNe enables us 
to address the question of the distance scale of these ob- 
jects. In Paper III (Napiwotzki & Schonberncr 1995| ) we 
proposed the use of the interstellar Na D lines for distance 
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Table 1. Model atmosphere fluxes (in 
10 8 ergcm _2 s^ 1 A~ 1 ) at A e g = 5454A for the calcu- 
lation of absolute magnitudes and distances. These 
values are calculated for a hydrogen-rich composition 
(niie/nu = 0.1) and logg = 6.5. However, the V band 
flux is quite insensitive to composition and log g 



in 10 erg cm" ^s" 1 !" 1 . We use the monochromatic ap- 
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determinations. These distances are free of assumptions 
about the nebula or the central star and are on average 
larger than the Shklovsky distances by a factor of 2.5 for 
our sample of old PNe. During the last years a number of 
trigonometric parallax measurements became availa ble for 
a sample of central star s of old PNe (Harris et al. 1997 , 
Pottasch & Acker |1998| , Gutierrez-Moreno et al. |1999| ). 
We will show that the CSPNe distances derived form our 
model atmosphere analysis, the NaD distances, and the 
trigonometric parallaxes are in agreement, but all three 
distances scales are much larger than those based on the 
Shklovsky method. 

2. Distance scales 

2.1. Spectroscopic distances 

Spectroscopic distances can be derived from the model 
atmosphere analyses described in Paper IV and offer an 
approach to the PN distance scale independent from the 
properties of the nebulae. After the stellar mass is esti- 
mated from a comparison with evolutionary models (cf. 
Paper IV) the distance can be calculated in a straight- 
forward way from effective temperature, gravity, and the 
dereddened apparent magnitude of the stars: 

1. the stellar radius is computed from surface gravity and 
mass, 

2. the surface flux in the V band is calculated from model 
atmospheres, with the appropriate parameters (thus 
bypassing the uncertain bolometric corrections), 

3. and finally the absolute magnitude My is derived and 
compared with the dereddend Vq magnitude of the 
star. 

4. This procedure can be condensed into one formula: 



d = o.ill jMlpL 10 oAv o 



(1) 



proximation of Heber et al. ( 1984j ) for Fy with an ef- 
fective wavelength A c ff = 5454 A, appropriate for hot 
stars. A short tabulation of model atmosphere fluxes 
is provided in Table |l|. 

This method is independent of any assumption about the 
PN, but obviously relies on the quality of the model at- 
mospheres used for the analysis and, to a lesser extent, on 
the evolutionary models used for the mass determination. 
For our objects the apparent magnitudes were taken from 
Table 4 of Paper III. If a high quality measurement (in- 
dicator "A" ) exists in the literature, it was adopted. Oth- 
erwise our spectrophotometric measurements were used. 
Since only very few reliable extinction measurements are 
available for our objects, we estimated the reddening from 
the galactic interstellar extinction model of Arenou et al. 
(1992). However, reddening is generally small and there- 
fore only of moderate importance for our conclusions. The 
resulting distances are listed as Canute in Table |^ together 
with the analysis results of Paper IV. Error limits were 
estimated from the analyses errors given in Paper IV, the 
estimated errors of the V measurements, and the redden- 
ing. The dominant error source is in most cases the gravity 
determination. 



2.2. Trigonometric distances 



Recently Harris et al. ( |1997| ) published trigonometric par- 
allaxes of 11 CSPNe. Eight of these stars belong to our 
sample (the value given for A 7 is actually only an up- 
per limit). Pottasch & Acker ( [1998D discussed HIP PAR- 
COS parallaxes of three CSPNe including PHL 932 which 
belongs to our sample. The resulting distance values are 
listed as du in Table ^ and are compared with the model 
atmosphere results in Fig. |]. The error limits were com- 
puted from the measurement errors of the parallaxes. 

As it is evident from Fig. [j] and the following figures the 
error distribution of distance determinations is highly non- 
Gaussian with many outliers, which can skew the common 
weighting techniques to erroneous values. The reason is 
that the usually (implicitly) adopted Gaussian probabil- 
ity distribution gives high weight to deviant points. There- 
fore we decided to minimize the absolute deviations, which 
corresponds to a double sided exponential probability dis- 
tribution and provides a more robust estimate (see dis- 
cussion in Press et al. 1992). Since the distance errors are 



where M* is the solar mass in M Q , g is the surface 
gravity in cm/s 2 and Fy is the model atmosphere flux 



highly asymmetric, we did the comparison with the par- 
allaxes, which have roughly symmetric error limits. Since 
we in all cases compare two measurements suffering from 
large uncertainties, we performed pro forma a linear re- 
gression with allowance for errors in both directions and 
with the intersection fixed at zero. The error ranges given 
below correspond to the la error of the mean. 

The measured trigonometric distances are always 
smaller than the NLTE distances (cf. Fig. [j]). The weighted 
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Table 2. Comparison of different distance determinations: from our NLTE analysis in Paper IV g?nlte, trigonometric 
distances dn of Harris et al. (1997), from the interstellar NaD line ^aD, and the Shklovsky distances efehki (the latter 
two are taken from Table 6 of Paper III). The kinematical ages calculated (ikin) from the spectroscopic distances are 
also provided. Effective temperature T c r , surface gravity g, and stellar mass M were taken from Table 1 of Paper IV. 
For further input data cf. Table 6 of Paper III 
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mean of the distance ratios r b s , computed as described 
above, amounts to 



r obs = = 1.55 + 0.29 



(2) 



Taken at face value this in concordance with the conclu- 



sions of Jacoby & Van de Steene ( |1995| ), who com pared 
trigonometric parallax measurements of Pier et al. (1993) 
with spectroscopic distances and concluded that the model 
atmosphere results overestimate the distances by ~40%. 
However, one must consider that trigonometric parallaxes 
are on the one hand the most direct method for distance 
determinations, but on the other hand are also subject to 
heavy biases, which are introduced by random errors of 
measurement and on average cause the trigonometric par- 
allaxes to be overestimated. This was discussed as early as 



1953 by Trumpler & Weaver. Lutz & Kelker (1973) were 



the first to provide a quantitative correction, which is a 
function of the relative measurement error of the parallax 
(Tn/n. The original Lutz-Kelker correction is valid in a 
strict sense only, if nothing else is known about the star, 
which is usually an unrealistic assumption. Different kinds 



if corrections have to be applied if one deals with a vol- 
ume limited sample, a magnitude limited sample, or if the 
"lumi nosity functi on" o f the parent population is known 
(Lutz 11983 ; Smith |1987| ). 



We took a different approach and performed a Monte 
Carlo simulation to derive biases caused by the selection of 
CSPNe for parallax measurements and the accompanying 
measurement errors. Our working hypothesis is that the 
spectroscopic distance scale is essentially correct, i.e. no 
systematic errors are present. The results of the Monte 
Carlo simulations are used to test if this hypothesis is 
compatible with observations. 

In our simulation CSPNe were randomly created ac- 
cording to a simple model of the Galactic stellar den- 
sity distribution and a theoretical post-AGB track. A 
spectroscopic analysis is simulated, i.e. random measure- 
ment errors are added. The star is selected for parallax 
measurement if the "spectroscopic distance" is below a 
threshold value <i max . A measurement error of the par- 
allax on is added to the true parallax and the resulting 
distance used for a comparison with the spectroscopic par- 
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Fig. 1. Distances computed from the results of the NLTE 
analysis compared to trigonometric distances. The error 
bars correspond to the measurement errors given in Ta- 
ble ||. The solid line indicates the average ratio of the 
NLTE distances and the trigonometric distances, the dot- 
ted line equality 



allax, if the "measured parallax" is larger than a threshold 
value II m i n . -/Vgampie stars are collected and the mean ratio 
rMc = dr ^' E is computed analogous to the procedure ap- 
plied for the observed data. This process is repeated iVtotai 
times and a probability distribution is evaluated. Sample 
probability distributions are given in Fig. |^ (these are dis- 
cussed below). A detailed description of the Monte Carlo 
simulations, the input parameters, and their standard val- 
ues is given in the Appendix. 

An investigation of the influence of the different input 



parameters on the results can be found in Table A. 2. The 



derived bias between spectroscopic and trigonometric dis- 
tances is most sensitive to the chosen threshold value d max 
for the pre-selection of stars, the parallax measurement 
error an and the threshold parallax II m j n . Our standard 
model assumes d,n ax = 1000 pc. Note that one star in the 



Harris et al. ( 1997 ) sample has a much larger value of 
the spectroscopic distance (A 74; ^nlte = 1700 pc). For 
on the mean of the error limits in Table ||, on = 0.95, 
was adopted. The threshold parallax was set to n m i n = 
1.33 (corresponds to A 74). Our standard sample size, 
Agampic — 8, corresponds to the observed sample and the 
simulation is carried out for iV tota i = 100000 samples. 

The first qualitative conclusion, which can be drawn 
from Table A. 2 and Fig. || is that the "measured" ratio 
tmc is always considerably larger than unity, although the 
underlying distance scales used in the Monte Carlo simu- 
lations are identical. Our standard model yields 



750pc 



= lOOOpc 




Fig. 2. Results of Monte Carlo simulations for several 
values of the maximum allowed "spectroscopic" distance 
rf max . The distributions are normalized to an integrated 
area of 1. For cosmetic reasons we used Atotai — 10 6 for 
this plot 



(we give the median value and the range which contains 
68.3% of all samples as la limit). Although our measured 
value r bs = 1.55 ± 0.29 is larger it is still in the la range 
of the standard simulation. If we increase d max to 1500 pc 
(keeping A 74 in mind) the result of the Monte Carlo simu- 
lation increases to tmc = 1-50 ±0.38 close to the measured 
value. The impact of modifying d max on the probability 
distribution is displayed in Fig. [2|. If we on the other hand 
remove CSPNe with c?nlte > 1000 pc from our sample in 
Table | (NGC 6720 and A 74) we derive r MC = 1.42 ±0.29 
close to the value predicted for d max = 1000 pc. 

A paradoxical effect results if the threshold parallax 
n m in is increased: the bias increases for increasing n m i n , 
if the spectroscopic threshold distance is left unchanged! 
E.g. for n min = 3.0 mas we derive tmc = 1-51 ± 0.35. 
The reason is that we introduce a strong selection effect 
for stars with much too large measured parallaxes. This 
effect is only overcome, if LT m i n is increased to even larger 
values (II m i n > 5 mas). 

In principle, the expected bias can be reduced by lower- 
ing c£ max and an. However, if we e.g. adopt d max = 750 pc 
and on = 0.7 mas our observed sample is reduced to 3 
stars! 

We conclude that the mean value of tmc derived from 
our observed sample (Table ||) is well within the range pre- 
dicted by Monte Carlo simulations for perfect agreement 
of both distance scales. Thus the trigonometric parallax 
measurements provide no evidence that the spectroscopic 
distances are in error, but confirm the spectroscopic dis- 
tance scale, instead! However, due to the statistical uncer- 
tainties we cannot provide a definitive proof of an agree- 
ment of better than about 20%. 



r MC = 1.32 ±0.25 



(3) 
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2.3. Distances based on interstellar lines and extinction 

Another method independent of assumptions about cen- 
tral star or nebula properties is based on interstellar 
lines. In Paper III we used the interstellar NaD lines 
at 5890/96 A for this purpose. The CSPNe are too hot 
to show any photospheric NaD lines and the PNe mate- 
rial is extremely dispersed. In their survey Dinerstein et 
al. ( 1995| ) detected circumstellar NaD lines in nine PNe. 
However, all investigated PNe were young objects, and 
if there is any neutral sodium at all in the circumstellar 
matter of the central star of old PNe, the column density 
would be very low and the resulting nebula contribution 
negligible. Naturally this method is restricted to stars at 
low galactic latitude. The distances derived in paper III 
are given as c?Na d in Table g. 

Our Na D distances in paper III were determined from 
the map of t he interstellar NaD line strength of Bin- 
nendijk ( 1952 ). One might ask, whether the distance scale 
adopted in this work is still valid. To our knowledge 
no more recent collection of interstellar NaD equivalent 
widths is available. The reason is that using the equiva- 
lent width is out of fashion, because resolved interstellar 
lines in high resolution spectra provide more information. 
However, for our faint CSPNe equivalent widths are still 
a useful tool. 



Binncndijk (1952) adopted spectroscopic distances of 
B stars, which were determined by the Yerkes group (Ram- 
sey 195C , Duke L951 , and unpublished distances provided 
by W.W. Morgan). It is save to assume that systematic 
differences between these authors are small. We performed 
a check of the Ramsey (1950) and Duke (1951) distances. 
For this purpose we selected a representative subsample 
of 20 stars from each collection with Stromgren uvbyfi 
measurements in the Hauck & Mermilliod ( 1998] ) cata- 
logue. Temperatures and gravities were calculated with 
an updated version (Napiwotzki & Lemke, in prep.) of 
the photometric calibration of Napiwotzki et al. ( 1993| ) . 
The new gravity calibration is based on accurate trigono- 
metric parallaxes measurements of B stars with HIP PAR- 
COS. Since the calibration doesn't cover O stars and su- 
pergiants, these stars were excluded from the comparison. 
Masses were derived by interpolation in the evolutionary 
tracks of Schaller et al. ( 1992| ) and distances d uv b y /3 com- 
puted as described in Section 2.1. Finally these distances 
were compared with the spectroscopic distances provided 
by Ramsey (|1950D and Duke fll95l|) and mean ratios were 
computed as described in Sect. 2.2. The results are 



a>uvby(3 
^Ramsey 
duvbyfi 
^Dukc 



0.92 
1.17 



(4) 
(5) 



or in other words the Ramsey distances are on average 8% 
longer, the Duke distances 14% shorter than the uvby(3 
distances. 
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Fig. 3. Distances computed from the results of the NLTE 
analysis compared to NaD distances. The solid line indi- 
cates the average distance ratio 



Although the scatter of the individual determinations 
is large we can conclude that altogether systematic differ- 
ences between the old Yerkes system and modern photo- 
metric distances, which are tied to accurate HIPPARCOS 
parallaxes, are small. 

Fig. |3| shows that the agreement of the NaD distances 
with the model atmosphere analysis is good. The scatter is 
easily explained by the spatially highly variable extinction 
in the galactic plane. The average ratio of both distance 
scales amounts to 



znlte 



iNaD 



1.18 ±0.22 



(6) 



The extinction distance method is an approach simi- 
lar to our NaD method. The interstellar reddening of the 
central star or the PN is measured and compared to a 
distance-reddening relation derived from stars in the an- 
gular vicinity of the PN. Since reddening and distances of 
many stars have to be determined, it is unfortunately a 
costly method and has been applied to few PNe only. Ex- 
tinction distances of two PNe of our sample were measured 



by Saurer et al. (g_995|). Results are 800 ± 300 pc for Sh2- 
188 and 800 ± 400 pc for HDW3. The Sh 2-188 distance 



is in good agreement with our estimate of IOOO^qq pc. 
Our HDW3 distance of 1500±^ pc is moderately larger 
than the extinction distance. Saurer noted that his upper 
limit maybe uncertain because the interstellar extinction 
reaches a plateau value at 1200 pc. However, the error bars 
overlap anyway. 
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2.4- Shklovsky distances 



The Shklovsky ( 1956 ) method allows the calculation of 
PN distances from the measurement of the recombination 
line H/3 and the angular diamete r. Th e distance can be 
computed from (see e.g. Pottasch 1984, p. 115): 



ffehkl = 218 "'" 



-0.18 



e 0.2 F 0| 9 0.6 



(7) 



with distance dshki in pc, flux F^p in ergcm~ 2 s _1 , the 
angular diameter in arc sec, the electron temperature 
t in 10 4 K, the filling factor e, and the mass of the ion- 
ized nebular matter M lon in Mq. The dependence of cfehki 
on t and e is only small, hence it suffices to use stan- 
dard values (e.g. t = 1 and e = 0.75). However, major 
problems are caused by the need to assume a mass M- lon 
of the (ionized) nebula. The classical Shklovsky method 
adopts a constant mass, usually M ion = 0.2Af Q . However, 
for the few PNe with independent distance determinations 
a trend of increasing mass with increasing radius is ob- 

1992]) . 



served, (see e.g. Maciel & Pottasch |198C; Cahn et al. ^_ 
Thus Macie l & Pottasch ( |1980D , Daub fll982j ), and Cahn 
et al. ( 1992] ) proposed modified Shklovsky methods em- 
ploying an empirical mass-radius relationship. However, 



the Daub (1982) and Cahn et al. (1992) mass-radius rela- 
tions adopted an upper limit for M- loni which means that 
for old, large PNe a constant mass of the order O.2M is 
applied, similar to the classical version. Recent approaches 



by Zhang ( |1995| ) and by Van de Steene & Zijlstra ( |1995| ) 
propose to use the correlation between the radio contin- 
uum brightness temperature and radius instead. However, 
this method needs an empirical calibration, too. In Ta- 
ble H we present Shklovsky distances dshki collected from 
the literature. Most values were taken from the catalogue 
of Cahn et al. ( 1992 ) . Please, consult Table 6 of paper III 
for detailed references. 

A comparison of the Shklovsky and the model atmo- 
sphere distances is shown in Fig |]. We distinguish be- 
tween "ordinary" white dwarf CSPNe of spectral type DA 
and DAO (filled symbols) and the hybrid/high luminos- 
ity objects A 43, NGC7094, Sh2-68, and DeHt2 (open 
circles) and the non post-AGB objects (open squares). 
The evolutionary history of both latter classes are likely 
very different from standard evolution and may cause very 
different PN properties. The comparison shows that al- 
most all Shklovsky distances of the ordinary white dwarf 
CSPNe are smaller than the model atmosphere distances. 
The average (weighted by the error limits of our analysis) 
amounts to 



^NLTE 



= 3.2 ±1.3 



hkl 



(8) 



where we have adopted an la error of all Shklovsky dis- 
tances of 30%. 

Since c?NLTE/^Shki is directly related to the ionized 
mass of the PN by the Shklovsky formula M- 101 
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Fig. 4. Distances computed from the results of the NLTE 
analysis compared to Shklovsky distances. Filled dia- 
monds: ordinary white dwarf CSPNe of spectral type DA 
and DAO, open circles: hybrid/high luminosity objects, 
open squares: non post-AGB objects. The dotted line in- 
dicates equality, the solid the average distance ratio for 
the ordinary white dwarf CSPNe 

interpret this as evidence that the ionized masses of the 
old PNe are much higher than the adopted typical mass of 
0.2Af Q . The distance scales would be in agreement if we 
increase the adopted PNe mass to 3.5M Q . However, the 
large scatter is a first indication that this is an oversimpli- 
fication of the real situation, as we will show below. Note 
in passing that the non post-AGB objects (open squares 
in Fig. H) have considerably lower values of cZnlte / <^Shki 
(mostly < 1) indicating masses lower than the adopted 
PN mass. 

3. Discussion 

We have now shown that the distances cZnlte derived from 
our NLTE model atmosphere analyses are (within today's 
error limits) in good agreement with the trigonometric 
parallaxes measured by Harris et al. (1997) and Pottasch 
& Acker (1998) after biases are taken into account, and 
with distances derived from the strength of the interstel- 
lar NaD lines. Both distance scales are model indepen- 
dent and thus demonstrate that the NLTE analysis are 
not subject to large systematic errors. On the other hand 



J2.5 



we 



Napiwotzki et al. (1999) have shown that state-of-the-art 
analyses of hot white dwarfs performed independently by 
different groups can yield surface gravities, which differ 
systematically by up to 0.1 dex which translates into a dis- 
tance error of 12%. Such errors of the model atmosphere 
distance scale would be compatible with the trigonometric 
and NaD distances of our stars. Much larger systematic 
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Fig. 5. The individual ratios of c?nlte and dshki as func- 
tion of the nebula radius. The meaning of the symbols is 
the same as in Fig. |J. The line indicates our fit (Eq. ^) 
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Fig. 6. Sames as Fig. [|but this time we compare the dis- 
tances from trigonometric and spectroscopic parallaxes du 
with dshki- The Ciardullo ct al. (1999) results are plotted 
as asterisks. For other symbols cf. Fig. [|. The line indi- 
cates the fit from Fig. [| 



errors can be excluded. The Shklovsky distance scale of 
old PNe is shown to be too short, most likely caused by 
an underestimate of the PNe masses. 

Further insight can be gained when plotting (Fig. [5]) 
the individual ratios c^nlte / <^shki as function of the PN 
radii (given in Table ||). A strong correlation is present. A 
simple (double logarithmic) linear fit results in 



log 



^NLTE 



0.388 + 0.935 log 



R 
pc 



(9) 



This indicates that the PN masses increase with increasing 
radii. The trigonometric measurements essentially confirm 
this result (Fig. ^|). The effect of the 35% shorter trigono- 
metric distances in (Fig. |]) is small, because moderate 



changes of the distances essentially produce a shift paral- 
lel to our fit line. 

How does our result compare with other tests of the 



Shklovsky distance scale? Ciardullo et al. (1999) deter- 
mined PN distances from main sequence companions of 
the CSPN and concluded that statistical distance deter- 
minations overestimate the PN distances. Stasihska et 
al. ( |1991fl and Pottasch & Zijlstra ( |1992| ) applied the 
Shklovsky method to PN in the bulge of our Galaxy. 
Stasihska et al. concluded that, besides considerable scat- 
ter, the mean Shklovsky distance reproduces the dis- 
tance of the galactic bulge well, while Pottasch & Zijlstra 
claimed that the Shklovsky distances are systematically 
too high. Do these results contradict our findings in Fig. 

Ciardullo et al. ( 1999 ) derived spectroscopic parallaxes 
of main sequence companions of 14 CSPN and compared 
their results (and 7 trigonometric parallaxes) with four 
different statistical distan ce scales based on th e pre scrip- 
tions of Cahn et al. ( |l 99$ Maciel & Pot tasch ( |1980D , va n 
der Steene & Zijlstra ( |1995| ) and Zhang ( |1995| ). These are 
variants of the Shklovsky method (Eq. ||) in which the ion- 
ized mass Mi on grows as a function of radius (Cahn et al. 
199$ Maciel & Pottasch |l98(i| ; Zhang |1995| ) or of the r adio 



brightness temperature of the PN. Cahn et al. ( |1992 ) as- 
sumed that the ionized mass of small, ionization bounded 
PNe grows with radius until an upper limit is reached 
for larger density bounded PN. This has the effect that 
these distance determinations for PNe in the sample dis- 



cussed by Ciardullo et al. (1999) are essentially classical 
Shklovsky distances with M- mn — 0.135M Q . We performed 



a small correction to the Cahn et al. (1992) distances to 
transform them to the com mon M j on — 0.2Af Q scale and 
added the Ciardullo et al. (|l998| ) data points to Fig. [| 
Two conclusions can be drawn: 

1. the Ciardullo et al. results agree very well with the cor- 
relation found for the spectroscopic distances in Fig. ^, 

2. the Ciardullo et al. distances are on average smaller 
than the Cahn et al. distances, indeed. The reason is 
that the PNe investigated by Ciardullo et al. are, on 
average, smaller than the PNe in our sample. 

PNe in the galactic bulge are very useful for several 
tests, because their distance is known from the distance of 



the bulge. Stasihska et al. (1991) and Pottasch & Zijlstra 
( 1992| ) used this for tests of the Shklovsky distance scale 
(both groups adopted a distance of 7.8 kpc). Stasihska et 
al. ( |1991 ) found reasonable agreement of Shklovsky dis- 
tance: their distance distribution peaked between 8 and 
9 kpc with 77% of all PNe in the range 8.5 ±3.5 kpc. Pot- 



tasch & Zijlstra (1992) on the other hand derived a peak 
value of 11.5 kpc and a 75% range from 7 to 17kpc. The 
main difference between both works is that Pottasch & 
Zijlstra used only fluxes and angular diameters measured 
from radio observations, while Stasihska et al. combined 
radio and optical flux data and preferred optical angular 
diameters. 
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Since one wants to exclude foreground objects, PNe 
with an angular diameter larger than 20" are excluded 
from bulge samples. This translates into a radius R = 
0.41 pc. Thus the old and large PNe of our sample (Ta- 
ble ^) are explicitly excluded from investigations of bulge 
PNe. Another strong selection effect against old PNe 
might already be at work, because of their low surface 
brightness combined with the large extinction. A lower 
limit of the angular diameter of w 1" is set by the need to 
resolve the PNe. This constrains the radii of bulge PNe, 
which could be used to test the Shklovsky method ap- 
proximately to the range 0.04 pc < R < 0.41 pc. From 
Fig. H we would predict that the Shklovsky distances of 
bulge PNe are moderately too larg e, in qualitative agree- 
ment with the Pottasch & Zijlstra ( 1992] ) result. However, 
one should keep in mind that the stellar population of the 
galactic bulge is quite different from the local population. 
Therefore one should be aware that the properties of bulge 
PNe might be different from local samples. 

The correlation given in equation ^ translates into a 
mass-radius relation 



1.87R 



2.3 



(10) 



with M in M and R in pc. Equation [L0| implies that the 
ionized masses of old PN are quite high. For an old PN 
with R = 2pc one derives Mi on = 9.3M . 

It is likely that during the PN evolution the com- 
plete material of the slow (fslOkm/s) AGB wind mate- 
rial is swept-up by the faster («20 . . . 40 k m/s) expanding 
planetary sh ell (Sc hmidt- Voigt & Koppen 1987 ; Marten & 
Schonberner 1991 ) and could finally be incorporated into 
the PN. However, 9M is much larger than the mass a 
typical CSPN precursor (M = 1. . .2M Q ) can loose dur- 
ing it's evolution. There are some points which should 
be considered when interpreting the mass-radius relation 
(Eq.0): 

1. the mass-radius relation would be modified if the filling 
factor e (cf. Eq. ||) changes during PN evolution, 

2. many PN of our sample (Table |J) have a large angu- 
lar diameter (several arcminutes are not uncommon) 
and the surface brightness is very low. Thus it is often 
a challenging observational task to measure emission 
line fluxes of these objects. As a consequence the ob- 
servational data provided in literature is often of low 
quality. We intend to determine reliable H{3 fluxes of 
old PN from Table ^ to improve this situation. 

3. Old PN show often signs of an interaction with the sur- 



rounding interstellar medium (Borkowski et al. 1990; 
Tweedy & Kwitter 1996 and references therein). The 
PNe shape is distorted and it is difficult to estimate 
the angular diameter O. 

On the other hand one may speculate that part of the 
ionized material is even swept-up interstellar matter. 
Note that an interstellar particle density of 1 cm~ 3 cor- 
responds to 0.05M Q /pc~ 3 . If the interstellar matter 




Fig. 7. Probability distribution for our standard model 
and for the possible future improvements discussed in the 
text 



corresponding to the nebular volume is incorporated 
into the nebula the mass of the R — 2pc PN is in- 
creased by 1.7 Mq. This might be an interesting topic 
for future investigations of the PNe interstellar matter 
interaction. However, more high quality observational 
material is necessary before one can start to attack this 
question. 

If the mass-radius relation in Eq. [l^ is confirmed one can 
use it in principle to produce a new variant of a statis- 
tical distance scale. However, the exponent (/3 = 2.3) is 
uncomfortable close to the value /3 = 2.5 at which the 
method degenerates and una ccept able large errors have 
to be expected (cf. e.g. Kwok |2000|, p. 177 f.). 



4. Conclusions 

We used the results of our model atmosphere analysis to 
compute distances of the central stars and showed that 
our model atmosphere distance scale is in good agreement 
with measured trigonometric parallaxes (after the effect 
introduced by biases has been taken into account via a 
Monte Carlo simulation) and distances derived from the 
interstellar NaD lines. All these three independent meth- 
ods of distance measurement indicate that the widely used 
"statistical" distance scales of the Shklovsky type are too 
short for old PNe. The most likely explanation is an un- 
derestimate of the nebula masses. 

Estimates of the PNe space density and birthrate, 
which are based on Shklovsky type distances, therefore 
give too large values. If a more realistic distance scale 
is applied, discrepancies between white dwarf and PNe 
birthrates are resolved. 

Due to statistical uncertainties and biases of the 
trigonometric distances we could test the spectroscopic 
distance scale only on the 20% level. Do we have to wait 
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Table A.l. 

simulation 



Standard parameters for the Monte Carlo 



Table A. 2. Influence of parameter variations on the re- 
sults of the Monte Carlo simulation 



Vmax 
rlmin 

on 



Sample sizes: 

iVtotal = 10 5 
^sample 8 

Sample selection: 

17.5 mag 
1000 pc 
1.33 mas 
0.95 mas 
0.75 mag 
Stellar evolution: 
Mcspn= 0.605 M© 
Tmin = 3000 yrs 
T max = 50000 yrs 
Galactic model: 
Ro = 4500 pc 
Hstar = 300 pc 
H dust = 140 pc 
Kdust = l.Omag/kpc 



for space missions like GAIA or SIM to improve this situa- 
tion? A test on the 10% level might be ach ieved with state- 
of-the-art-techniques. Harris et al. ( 1997 ) announced that 
parallax measurement accuracies of 0.5 mas are within 
reach, but this alone is not enough (cf. Table A. 2 ). How- 
ever, if a sample with c?nlte < 750 pc is produced and 
the sample size is doubled to N samp \ c = 16 the bias rue is 
reduced to 4% and the scatter to 10% (the probability dis- 
tribution is compared with the present situation in Fig. Q). 
Thus, if we work on both fields, measuring of more accu- 
rate parallaxes and analyzing and selecting more candi- 
dates from spectroscopic investigations, this goal is within 
reach. 

Acknowledgements. The author thanks Joachim Koppen for 
inspiring discussions and Detlef Schonberner, Uli Heber, and 
Klaus Werner for useful comments on previous drafts of this 
paper. 

Appendix A: Monte Carlo simulation of the local 
PN population 

Our Monte Carlo simulation of the local PN distribution pro- 
ceeds in three steps: 

1. a star is created at a random position based on a simple 
3D model of the galactic disk, 

2. a post-AGB age is randomly chosen and the absolute mag- 
nitude of the star is computed from a theoretical evolution- 
ary track, 

3. a scheme which describes the selection of CSPNe for par- 
allax measurements and the effects of measurement errors 
is applied. 

The par ameters of our standard model are provided in Ta- 
ble A.l and the results of a parameter study which exp lores 



parameter 


values 


^NLTE/rfn 


std. model 




1.32±0.25 


N sample 


16 


1.30±0.19 




32 


1.30±0.13 




64 


1.30±0.09 


Vmax 


16.5 mag 


1.24±0.25 




18.5 mag 


1.34±0.25 


^max 


750 pc 


1.14±0.20 




1500 pc 


1.50±0.38 




3000 pc 


1.64±0.51 


n m j n 


1.5 mas 


1.37±0.26 




2.0 mas 


1.51±0.30 




3.0 mas 


1.51±0.35 




4.0 mas 


1.32±0.30 




5.0 mas 


1.22±0.24 




7.0 mas 


1.16±0.21 




10.0 mas 


1.14±0.20 


cm 


0.1 mas 


1.06±0.16 




0.2 mas 


1.07±0.16 




0.3 mas 


1.10±0.16 




0.4 mas 


1.13±0.17 




0.5 mas 


1.16±0.17 




0.7 mas 


1.23±0.20 




1.2 mas 


1.38±0.31 


o~m v 


0.5 mag 


1.26±0.23 




1.0 mag 


1.30±0.28 




1000 yrs 


1.31±0.26 




5000 yrs 


1.31±0.26 


Tmax 


25000 yrs 


1.33±0.26 




100000 yrs 


1.29±0.26 




150 pc 


1.28±0.26 




600 pc 


1.34±0.26 


Hdust 


60 pc 


1.32±0.25 




210 pc 


1.31±0.25 


^dust 


0.50 mag/kpc 


1.32±0.26 




1.50 mag/kpc 


1.30±0.26 



Our simple description of the galact ic disk is based on the 
Galaxy model of Bienayme et al. (1987). An exponential den- 
sity law with a scale height // star of 300 pc is adopted for the 
CSPN. This corresponds to a stellar population with an age 
of 3. . .5 • 10 9 yrs. We included extinction by a dust compo- 
nent with a scale height Hdust = 140 pc, which is the value 
appropriate for the interstellar matter. A dust opacity at the 
position of the sun ftdust ~ 1 mag/kpc was chosen. The stel- 
lar and dust density decreases exponentially with the distance 
from the galactic center and a scale length Ro = 4.5 kpc. For 
the distance of the sun from the galactic center we adopted the 
stand ard value of 8.5 kpc. Let us note that Koppen & Vergely 
(1998) could successfully reproduce the extinction properties 



the effects of varying these parameters are listed in Table A. 2 



of galactic bulge PNe with our parameter values. 

Our parameter study in Table A. 2 shows that the influence 
of particular values of the parameters of our "galactic model" 
are very small. Due to the exponential decrease of stellar den- 
sity with height above the galactic plane the number of stars 
within a sphere with a given radius R increases less then R 3 . 
Since this lower the numbers of far away stars, a lower bias 
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is expected for lower values of the scale height. Extinction in- 
troduces another selection against far aw ay st ars (through the 
limiting magnitude). However, as Table A. 2 proofs the effect 
of varying these parameters within reasonable limits is quite 
small. 

For a given post-AGB age the absolute magnitudes of 
CSPNe can be computed from theoretical post-AGB tracks. 
In principle a mass distribution for the CSPNe should be 
used. However, since the mass distributions of CSPN and white 
dwarfs are quite narrow, we considered it sufficient to use only 
one theoretical track for simplicity. Spectroscopi c stu dies of 
CSPNe (Pap er IV) , white dwarfs (Bergeron et al. |l992j , Napi- 
wotzki et al. 1999) and an investigation of P Ne bas ed on a 



distant independent method by Stasiriska et al. ( 1997 ) derived 
peak masses in the range 0.55. . .0.60Mp, . Thu s we selected the 
O.6O5M post-AGB track of Blocker (|l995|) for our simula- 



tion. The post-AGB age was varied within a given interval 
Tmin- • -Tlnax which approximately reproduces the My distri- 
bution of our CSPNe. However the simulation resul ts are quite 
insensitive to their particular values (cf. Table A.2). 

After "producing" a CSPN we had to simulate the selection 
for parallax measurement and the distance determination with 
their measurement errors. The following scheme was adopted: 

1. we simulated a "spectroscopic analysis" with a random 
measurement error. We adopted a standard value for the 
scatter of the absolute magnitude determination of om v = 
0.75 (Gaussian error distribution). This corresponds to the 
mean of the errors given in Table |^. 

2. if the "spectroscopic analysis" indicates that the distance 
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of the CSPN is helow a given maximum value d m?i .. the star 



ph/9908181 



astro- 



was selected for "parallax measurement". 

3. the parallax is "measured" with an Gaussian error distri- 
bution with an- Our standard value of an is 0.95 mas. That 
is the mean error of the parallax measurements provided in 
Table ^. A parallax measurement is used for distance deter- 
mination if the value is larger than a lower limit n m j n . We 
chose n m i n = 1.33 mas, corresponding to the value mea- 
sured for A 74. 

4. this process was repeated until iV aam pie distance determi- 
nations were performed. The standard value of -/V aamp i c is 
8, the number of CSPNe with parallax measurements in 
Table |^. The mean value of g(nlte / dn was computed with 
the same method as used for our observed sample (cf. Sec- 
tion 2.2.). 

5. to derive a proper statistic this simulation was repea ted 
-Wtotai = 10 s times. The best values given in Table A.2 and 
in Section 2.3 are the median values, la values were com- 
puted from the condition that they should include 68.3% of 
the simulations. Results from repeated runs yielded results, 
which didn't deviate by more than 1%. 

The effect of the most important input parameters, the thresh- 
old value d m ax for the pre-selection of stars, the parallax mea- 
surement error an and the threshold parallax n m j n , is discussed 
in Section 2.2. 
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